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We study the dynamics of a beam of fermions diffracted off a density grating formed by fermionic 
atoms in the limit of a large grating. An exact description of the system in terms of particle-hole 
operators is developed. We use a combination of analytical and numerical methods to quantitatively 
explore the Raman-Nath and the Bragg regimes of diffraction. We discuss the limits in diffraction 
efficiency resulting from the dephasing of the grating due to the distribution of energy states occupied 
by the fermions. We propose several methods to overcome these limits, including the novel technique 
of "atom echoes". 
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I. INTRODUCTION 

The availability of quantum-degenerate bosonic atomic 
systems has recently allowed the extension of atom optics 
to the nonlinear [n| and the quantum regimes. Matter- 
wave four- wave mixing , coherent matter- wave ampli- 
fication 1^ and superradiance , the generation of dark 
1^ and bright ^ atomic solitons and of correlated atomic 
pairs have been demonstrated, and so has the matter- 
wave analog of second-harmonic generation, the creation 
of a molecular condensate component || . 

In contrast, the development of the nonlinear atom op- 
tics of fermionic atomic systems is not nearly as far along. 
While it has been shown theoretically |ll| that the 
four-wave mixing of fermionic matter waves is possible 
in principle, these predictions have not been verified ex- 
perimentally so far. Still, the recent achievement of tem- 
peratures as low as 0.2Tp, where Tp is the Fermi tem- 
perature, for the fermions and ^Li [l3[ |l| 
is encouraging, and it is hoped that first experiments on 
fermionic nonlinear atom optics will take place in the near 
future. In addition to the fundamental goal of achiev- 
ing a BCS phase transition into pairing and superfluid- 
ity research along these lines is also motivated by 
recent results that hint at the possibility to lower the 
phase noise in intcrfcromctric measurements below the 
bosonic standard quantum limit by using instead degen- 
erate fermionic beams |p7| . 

The first theoretical discussions of fermionic nonlin- 
ear atom optics were presented in Refs. and [pd]| , 
which treated the case of a single 'test' particle scatter- 
ing off a periodic density grating formed by a degenerate 
Fermi gas. They showed that for an appropriately pre- 
pared grating, the fermionic system can indeed undergo 
four-wave mixing. In contrast to the standard interpreta- 
tion in terms of "bosonic amplification" , which clearly is 
not applicable to fermions, this effect was interpreted in 
terms of the constructive quantum interference between 
different paths leading to the same final state. 

One important aspect of the fermionic case is that, in 
contrast to bosons, considerable care must be taken in 
combining two matter waves to form a "grating" , so that 
their interaction with a third wave can produce a signifi- 



cant four- wave mixing signal. Consider, as we shall do in 
this paper, counterpropagating matter waves separated 
in momentum by q. In the case of bosons, two obvious 
possibilities correspond to the states 
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the 6 J, being usual bosonic creation operators and |0) 
the atomic vacuum. The fist case describes two counter- 
propagating beams of N/2 atoms each and of momenta 
±g/2, while the second state corresponds to a density 
grating obtained by identically preparing N atoms in a 
coherent superposition of states of momenta ±g/2. It is 
known from the study of atomic diffraction by optical 
fields that these two states lead to different diffraction 
patterns, because the first one contains "which way" in- 
formation while the second doesn't This difference 
becomes however insignificant for large gratings. 

The situation for fermions is more complicated, since 
the Pauli exclusion principle precludes one from placing 
more than one atom per state. One needs instead to 
consider multimode atomic beams, centered around the 
mean momenta ±q/2. In this case the states \tpb,i) and 
|V'fc,2) are replaced by 
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where a[._(_q/2 fermionic creation operators for atoms 
of momenta in the vicinity of the total number of 

atoms involved being indicated in the appropriate prod- 
ucts. From Refs. |l^, we know that it is the 
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FIG. 1: Spin-down polarized fermionic beam diffracted by a 
second species of spin-up polarized fermions forming a matter- 
wave grating. 



quantum coherence apparent in matter-wave states of 
the form |'0/,2) that is responsible for fermionic four-wave 
mixing. In order for the required quantum interference 
to occur, it is essential that every atom be in a coherent 
superposition of momentum states centered around g/2 
and ~q/2. 

So far, our discussion has ignored the time dependence 
of the grating. But since the atoms forming a fermionic 
grating all have slightly different kinetic energies, their 
free evolution results in a dephasing that is expected 
to eventually lead to the disappearance of the four-wave 
mixing signal. Although the importance of this dephas- 
ing was pointed out in Rcf . [p^ , no quantitative discus- 
sion of its effect has been presented so far. The present 
paper addresses this problem quantitatively by a com- 
bined analytical and numerical study of the diffraction 
of a beam of fermionic atoms off a large fermionic grat- 
ing. We fully include the dynamics of the atomic beam, 
but neglect its back-action on the grating dynamics, con- 
sidering only its free evolution. This is the matter-wave 
analog of the undepleted pump approximation in nonlin- 
ear optics. 

Section II introduces our model and describes its dy- 
namics in terms of equations of motion for particle-hole 
operators. The effect of the grating dynamics is charac- 
terized in terms of a dephasing time r^;, whose impact is 
then illustrated in the simple case of Raman-Nath diffrac- 
tion. The Bragg regime is analyzed in section III using a 
combination of an analytical discussion of a simple two- 
mode system and full numerical simulations. We deter- 
mine the characteristic time scales governing this process, 
and conclude that four-wave mixing in degenerate Fermi 
gases should barely be observable. Noting that the de- 
phasing of the grating is in principle reversible, we turn 
in Section IV to a discussion of possible ways to achieve 
such a reversal, based on analogies with the photon echo 
techniques of quantum optics. Since the physical origin 
of the dephasing is the difference in kinetic energies for 



atoms with different momenta, one elegant approach ex- 
ploits the interaction of the grating with a periodic exter- 
nal potential to achieve a negative effective mass of the 
atoms. Such methods pave the way for the observation of 
phenomena that seem out of experimental reach at first 
glance. Finally, Section V is a summary and outlook. 



II. MODEL 

A typical atomic four-wave mixing experiment involves 
three input wave packets that interact, typically via s- 
wave collisions, to produce a fourth packet in a new mo- 
mentum state. Physically, this can be interpreted as the 
result of the phase-matched scattering of one of the wave 
packets off the matter-wave grating produced by the su- 
perposition of the other two, as discussed in the intro- 
duction. 

As is readily apparent from Eq. (^, though, a grating 
formed by N fermions occupies 2A^ momentum states, 
with TV states around the momentum q/2 and N states 
near —q/2 in our case. As a result of the free evolution 
associated with the spread in kinetic energies of these 
atoms, the density grating undergoes a dephasing, similar 
to Doppler broadening in a thermal vapor. The goal of 
this paper is to investigate the effects of this dephasing 
and to determine to what extent it can be reduced, or 
even eliminated. 

We consider the situation of a spin-polarized fermionic 
beam diffracted by an atomic grating of spin-polarized 
fermions that are either of a different species or in a 
different spin state, see Fig. ^ the beam and grating 
atoms interacting coherently via elastic s-wave scatter- 
ing. We neglect p-wave scattering, whose cross-section 
at zero temperature is orders of magnitude smaller than 
that of s-wave scattering. 

The diffraction of a fermionic atomic beam by a pe- 
riodic potential was previously treated for the case of 
free fermions (or equivalently a one-dimensional bosonic 
Tonks gas) |^ and of Cooper pairs j2^. These treat- 
ments are restricted to the Raman-Nath regime, where 
the kinetic energies of the diffracted atoms are negligible 
compared to the potential energy of the grating. The 
scattering of two fermions by an optical standing wave 
was also considered in Ref. |22[. Also worth mentioning 
in the present context is Ref. [p3| , which investigates the 
nonlinear mixing of two bosonic and one fermionic wave. 
In these cases Pauli blocking has a significant effect on the 
diffraction pattern, but as the gratings are non- fermionic, 
the dephasing, which is central to our work, is not ad- 
dressed. 

We restrict the quantized description of the system to 
the dimension along the axis of the grating only, the 
"transverse" direction. The dynamics of the system in 
the direction perpendicular to it is treated classically, 
and parameterized by the time the atoms in the beam 
interact with the grating. Note that this approximation 
neglects both the normal reflection and the Bragg reflec- 
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tion of the incident beam off the grating ||2^ . 

This model is described by the second-quantized 
Hamiltonian 



H = 



fc 



(5) 



,(5) 



f. — hk^/2mg (i^^^"* = hk'^/2mb) is the kinetic 



where w 

energy of a grating (beam) atom of mass rUg {mt) and 
transverse momentum fc, 



Uo = 2TTha/{LAfi), 



(6) 



and k + q/2, with —kF,g < fc < kp^g, where fcp.g = T^Ng/L 
is the Fermi momentum (in one dimension) and Ng is the 
number of atoms forming the grating. Such a grating 
can be prepared with two-photon Bragg pulses of band- 
width larger than the Fermi frequency hkp g/ (2mg) in a 
two-step process starting from a stationary homogeneous 
gas, see e.g. Ref. [|4|. In that scheme, a 7r-pulse impart- 
ing a transverse momentum —hq/2 to each atom in the 
static cloud is followed by a 7r/2 Bragg pulse that creates 
the desired superposition state. Thus the initial state of 
beam-grating system is 



where 



(9) 



fi — TngTrifj/ {vrLg + mh) is the reduced mass, a is the s- 
wave scattering length characterizing the ultracold col- 
lisions between grating and beam atoms, A is an effec- 
tive transverse cross-sectional area, and L is the quanti- 
zation length. The annihilation and creation operators 
ak.s and at satisfy the fermionic anticommutation re- 



lations [afcs,afc/s/] , = and 



Skk'Ss 



the 



spin components s = {Tii} corresponding to the differ- 
ent spin states or atomic species of the grating and the 
beam, respectively. 

Since the elementary process underlying matter-wave 
diffraction consists in annihilating an atom of momentum 
k and creating an atom with momentum k + q, it is useful 
to describe the dynamics of the atomic beam in terms of 
"particle-hole" operators 
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where the subscripts "b" and "g" label beam and grating 
atoms, respectively. The Heisenberg equations of motion 
for these operators are easily found as 
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A similar equation holds for 



the grating particle-hole operators, with the substitution 
g. 

We consider a beam of Ni, atoms propagating toward 
the grating with some central transverse momentum q. 
As discussed in the introduction the grating is taken to be 
in a superposition of states of transverse momenta k—q/2 
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The corresponding initial expectation values of the 
particle-hole operators are 



pS'(o)|*^''^) = '5,,oe(fcF,fc-|fc-g|), 
pS(0) 
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In the following we assume that the matter-wave grat- 
ing is sufficiently large that one can safely ignore the 
back-action of the atomic beam on its dynamics. Math- 
ematically, this amounts to neglecting the term propor- 
tional to Uq in the grating version of Eq. (||). In this 
limit, the grating simply undergoes a free evolution fully 
governed by its transverse kinetic energy. 

Assuming, consistently with the assumption of a large 
grating, that the expectation value of products of beam 
and grating particle-hole operators can be factorized as 



(12) 



making explicit use of the grating initial condition (|l]^ 
and introducing the slowly varying particle-hole opera- 
tors 

pS'W-/5S'Wcxp[-*4'J'^], 

and similarly for the grating, yields then the beam 
particle-hole Heisenberg equations of motion 
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FIG. 2: Time evolution of the transverse mode populations 
of the diffracted beam, (rife) = {Pk,o), in the Raman-Nath 
regime. The initial momentum spread of the beam is kp.b ~ 
3q. The effects of the dephasing grating become apparent. 
Time in units of the inverse mean-field interaction strength 
L^oA^g and — 15. 



where 



(^) = V E ©(^^.s - + 9/21) exp[-4^|^i], (14) 



reduces in the continuum Umit, — > L/2t: J dk, to 
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with sinc(.T) = 
phasing time 



sin(a;)/x. For times larger than the de- 



hqkp.g ' 
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we have that g{t) ~ 0, that is, the diffraction of the beam 
comes to a stop. This consequence of the dephasing be- 
tween different momentum states of the atomic grating 
resulting from their free evolution, represents an essen- 
tial difference between fermionic and bosonic four-wave 
mixing: the free evolution of a bosonic grating optically 
prepared from a Bose-Einstein condensate at zero tem- 
perature gives rise to a phase factor that can easily be 
transformed away. For fermions, though, it results in a 
dephasing and the effective shutting off of its interaction 
with the atomic beam. 

This impact of the dephasing is particularly easy to 
see in the Raman-Nath regime of atomic diffraction. We 
recall that this is the regime where the kinetic energy of 
the atoms in the beam is small compared to the poten- 
tial energy of the density grating. However, we retain 
the kinetic energies of the grating atoms, since they are 
responsible for the dephasing. As such, this model corre- 
sponds to the case where the atoms in the incident beam 
are much heavier than the atoms in the grating. 



Approximating the exponential functions by unity in 
Eq. (^3|) consistently with the Raman-Nath approxima- 
tion gives 
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The solution of these operator equations for g{t) = con- 
stant are known j20| in terms of Bessel functions of the 
first kind of order s as 
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Due to conservation of momentum, only states of trans- 
verse momenta separated by integer multiples of q are 
dynamically coupled. For a beam of initial momentum 
spread kp^b < 9/2, this implies that a linearly increas- 
ing sequence of replica of the initial distribution will 
be generated in time. Things are more complicated if 
kF,b > 9/2, since Pauh blocking leads to a situation where 
only those transverse modes near the edge of the distribu- 
tion can initially be diffracted, creating holes into which 
modes deeper into the initial distribution can then sub- 
sequently be coupled. This results in an effective broad- 
ening of the initial distribution in time, as illustrated in 
the early stages of Fig. ^. 

Since g{t) is actually not a constant, the solution ( |l7| ) is 
only valid for times t <^Td- In general, it is necessary to 
solve the expectation value of the coupled Raman-Nath 
equations numerically. The effect of the grating diffusion, 
and the resulting vanishing of g(t) for long enough times, 
are readily apparent in Fig. ^, which shows that the beam 
diffraction effectively ceases for t > Td- This is seen even 
more clearly in Fig. ^, which compares the number of 
significantly occupied transverse modes with and without 
dephasing of the grating. 



III. BRAGG REGIME 

We now turn to the analysis of the Bragg regime of 
fermionic atomic diffraction. We recall that for single 
atoms and a grating of period in/q, this regime corre- 
sponds to the situation where only two degenerate modes, 
of transverse momenta ±g/2, are coupled. The probabili- 
ties of occupation of these two states oscillate periodically 
in time, the so-called Pendellosung |2^. Bragg scattering 
is widely used in atomic beam splitters, both for single 
atoms and for Bose-Einstein condensates ||2^ . 

Fermi statistics complicate the situation in two ways: 
First, as we have seen, the distribution of momentum 
states in the grating around ±g/2 leads to a dephasing 
resulting in the time-dependent coupling constant g(t). 
In addition, only those fermions with initial transverse 
momentum sufficiently close to q/2 fulfill the Bragg con- 
dition. Hence, it is expected that Bragg diffraction will 



FIG. 3: Largest transverse momentum state with a signifi- 
cant population (> .03) as a function of time in the presence 
(circles) and absence (crosses) of dephasing. Time is in units 
of the inverse mean-field frequency UoNg and Td — 15. 



FIG. 4: Time evolution of the transverse mode populations 
of a fermionic ®Li-beam scattering off a *°K-grating in the 
Bragg regime. Time is in units of the inverse mean-field 
frequency NgUo- In this example Td ~ 4.2/ (Ngllo) and 
LOq'''' / NgUo = 5. The full momentum widths of the grating 
and of the beam are l.Oq and 0.25q, respectively. 



burn a spectral hole in the initial Fermi sea of the beam, 
leaving atoms away from that hole practically untouched. 

We can gain some insight into the effects of the dephas- 
ing by considering a simple two-state model that retains 
only the coupling between the transverse modes 
neglecting their coupling to higher-order modes. This 
approximation is justified if the phase-matching condi- 
tion is strongly violated for these modes, in other words, 
if the difference in kinetic energy q much larger 
than the mean-field energy hUoNg responsible for mode 
coupling, i.e.. 
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q > Anapg , 



(18) 



Pg being the atomic density in the grating. 

This situation is described by the closed set of operator 
equations 



d_ .(b) 



±ig{t) 



4b) 

,ib) 



«/2 



^-1/2 



(19) 



where g{t) is given by Eq. ([l5|). These equations are 
reminiscent of the optical Bloch equations that describe 
the dynamics of a two-level atom driven by a resonant 
single-mode field of time-dependent amplitude, with the 
operators n corresponding to the populations of the two 
levels and the particle- hole operator q atomic 
polarization. The evolution of their expectation values 



can readily be found for the initial conditions 
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where 



Pit) 



dt'g{t') = -NgUo / dt'sinc 



hkp^gq 



Similarly to the situation that we encountered in the 
Raman-Nath regime, the dephasing of the grating results 
in the diffraction effectively coming to an end after a time 
of the order of t^. 

Fig. ^ summarizes the results of the numerical inte- 
gration of the full equations of motion (^) in the Bragg 
regime for a ^Li beam scattering off a grating formed by 
atoms. It shows the population of the transverse 
modes of the atomic beam as a function of time for the 
case where the initial full transverse momentum width 
of the grating and beam are l.Oq and 0.25(7, respectively. 
The hole burned in the initial Fermi sea is clearly vis- 
ible, as is the impact of the dephasing of the grating. 
The first maximum in the Pendellosung oscillations is 
strongly pronounced for modes near the phase matched 
mode q/2, but further maxima cannot be observed due 
the dephasing of the fermionic grating. 

One can characterize the "efficiency" rj of the four- wave 
mixing process by the first maximum of the quantity 



(22) 



which measures the maximum fraction of atoms whose 
transverse momentum has been flipped about fc = 0. It 
is plotted in Fig. [q as a function of the momentum spread 
of the incident Li beam for three widths of the "^''K 
grating. The behavior of 77 is as expected, since a wider 
incident beam contains more atoms that do not fulfill the 
Bragg condition, and a larger Fermi sea of the grating 
means a shorter dephasing time. 
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2 ^/q 

FIG. 5: as a function of beam size, kp^i,, for different sized 
gratings, with ulj''' /Ngllo = 5. The efficiency rj decreases with 
separately increasing size of the Fermi seas of beam and grat- 
ing. Nevertheless for reasonably small momentum spreads a 
decent efficiency is obtained. 
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FIG. 6: Effective masses of ''"K of the first two bands at the 
center of the first Brillouin zone as a function of Vo, in units 
of the recoil energy iJ,, — h^k^ /2mg for an optical wavelength 
of 852 nm. 



Experiments are likely to be carried out in a three- 
dimensional geometry, in which case the time-dependent 
coupling constant (|lj) becomes 

53Z)(t) = 3UoNg (sine— - cos — ) rj/t^, (23) 

where Td is still given by Eq. (p^. It is easily seen that 
gsoit) decays in a manner similar to g{t), so that the 
lifetime of the grating is essentially the same as in one- 
dimension, provided that kF,g remains the same. Note 
however that these two times scale differently with the 
density of the grating, due to the different mode densities 
in one and three dimensions. 

In order to assess the experimental feasibility of a 
fermionic Bragg diffraction experiment it is necessary 
to compare the dephasing time Td to the time tb that 
it would take to observe a full Pendellosung oscillation 
in the absence of dephasing. With the approximation 
f3{t) « NgUot/2 we have readily that tb ~ ir/NgUo- 
Hence, the observation of Pendellosung oscillations re- 
quires 

where we have used the three-dimensional relationship 
between the density pg and kp.g, kp.g = (GTr^pg)^/^. 

Consider for example the diffraction of fermionic potas- 
sium ^"^K off a spin-polarized ^''K grating of peak number 
density pg — 10^^ cm^'^ and a scattering length of 85ao. 
Assuming a grating period of 500 nm Td ~ O.Itb, so that 
Pendellosung oscillations will not be observable in that 
case. The next section discusses possible ways around 
this difficulty. 



IV. ATOM ECHOES 

The condition (^J) shows that one possible way to in- 
crease the ratio of Td to tb is to increase the grating 
density pg. This, however, is likely to be impractical. An- 
other solution is to increase the scattering length a via a 
Feshbach resonance. Noting that rfigj p = (mg + nih)/mb, 
we can also gain an order of magnitude or so in the ratio 
( p^ ) by diffracting light fermions such as ^ Li off a potas- 
sium grating. Ketterle et. al Q also proposed a direct 
way to increase the decay time in three dimensions by 
reducing the effective Fermi wave vector in the direction 
of the grating by increasing the size of the Fermi sea in 
the other directions. 

An alternative method that allows one to completely 
reverse the effects of the dephasing is based on the obser- 
vation that this dephasing is coherent, simply resulting 
from the different kinetic energies of the atoms forming 
the grating. This is reminiscent of Doppler broadening 
in laser spectroscopy, which leads to a dephasing of the 
atomic polarization. In that case, the dephasing can eas- 
ily be reversed by photon echo techniques ||2^. Things 
are slightly more complicated here, because the phase of 
the individual atoms is just given by their kinetic en- 
ergy h^k"^ /2mg. The only way to reverse the sign of this 
phase is to reverse the mass of the atoms. This can be 
achieved by placing the atoms forming the grating in an 
additional periodic potential that gives them a negative 
effective mass. 

Ignoring two-body collisions as well as whatever exter- 
nal fields are required to prepare the initial grating ( |l0| ) 
for now, the evolution of the grating atoms in the pres- 
ence of an optical lattice is governed by the Hamiltonian 

H^^ + ^cos{2kLx), (25) 

where Vb is the lattice depth and fcz, the wave vector of 
the laser creating the periodic potential. The band struc- 



ture associated with this potential is readily obtained nu- 
merically, yielding at the center of the first Brillouin zone 
the effective mass 
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(26) 



where n is the band index and k is the quasi-momentum. 
Figure |^ plots this effective mass of '*°K for the first two 
energy bands. As is well known, is alternatively 

positive and negative as the band index is increased. We 
remark that since Uo is a pseudo-potential adjusted so 
that the scattering problem associated with two-body 
collisions asymptotically produces the correct coUisional 
cross-section, the appearance of an effective mass does 
not change the reduced mass in Eq. (||). Hence Eq. 
(U) becomes 
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From Fig. ^, we see that the effective mass can be orders 
of magnitude larger than nig , leading to much easier con- 
ditions to satisfy to observe fermionic Bragg diffraction. 

Even more interesting is the use of the negative effec- 
tive mass associated with the second energy band, since 
it permits a full reversal of the dephasing. The excita- 
tion of the grating atoms to that band can be achieved 
e.g. via a two-photon Raman with co-propagating laser 
beams, or by modulating the phase of one of the lat- 
tice beams. Assuming for simplicity that the effective 
masses corresponding to the n — 1 and n = 2 bands are 
equal in magnitude, it is easily seen that if the Raman 
pulse is applied at time ti, then g{t2 = 2ti) — g{0), and 
the dephasing has been completely reversed. Clearly, a 
sequence of pulses alternatively transferring the grating 
atoms between the two energy bands will keep alternat- 
ing the sign of rrieff, thereby continuously compensating 
the Doppler broadening. In case the magnitude of curva- 
tures of the two bands are different, the rephasing time 
needs to be appropriately adjusted. 

Fig. 1^ shows the result of this atom echo technique. 
The second peak of the Pendelldsung is now clearly visi- 
ble, in contrast to the situation in Fig. ^. For better com- 
parability to the case of uncontrolled dephasing we con- 
sidered the same magnitudes for effective and real mass. 
We have transferred the atoms to the first band at scaled 
time t = 3 and back to the lowest band at t = 9. The re- 
sults of the dephasing of the grating start again later on, 
but could be eliminated by transferring the atoms back 
to the first energy band. 



V. CONCLUSION 

In this paper, we have investigated the behavior of 
purely fermionic scattering, i.e. a beam of fermions 
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FIG. 7: Comparison between the rephasing of the Pen- 
dellosung oscillations by the atom echoes technique (upper 
plots) to the case where the dephasing has been arbitrarily 
turned off (lower plots). In the plots white corresponds to 
unit probability and black indicates zero probability. All pa- 
rameters as in Fig. ^. 



diffracted by a fermionic density grating, in the limit of 
a large grating. We explored different regimes which are 
determined by the relative magnitude of mean field en- 
ergy and kinetic energy. Fermi statistic shows its distinct 
signatures in the limited lifetime of the grating and the 
inability to satisfy the Bragg condition for large input 
beams. We proposed novel ways to increase the dephas- 
ing time of the fermionic grating which is the most se- 
vere limitation for fermionic four-wave mixing. Future 
work will focus on the complex dynamics arising from 
the back-action of the scattered atoms on the grating. 
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